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MULTIADAPTIVE GALERKIN METHODS FOR ODES III: 
A PRIORI ERROR ESTIMATES* 
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Abstract. The multiadaptive continuous/discontinuous Galerkin methods mcG{g) and mdG(q) 
for the numerical solution of initial value problems for ordinary differential equations are based on 
piecewise polynomial approximation of degree q on partitions in time with time steps which may vary 
for different components of the computed solution. In this paper, we prove general order a priori 
error estimates for the mcG(g) and mdG(g) methods. To prove the error estimates, we represent the 
error in terms of a discrete dual solution and the residual of an interpolant of the exact solution. The 
estimates then follow from interpolation estimates, together with stability estimates for the discrete 
dual solution. 
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1. Introduction. This is part 3 in a sequence of papers |321[S3] on multiadap- 
tive Galerkin methods, mcG{q) and mdG(<7), for approximate (numerical) solution of 
ODEs of the form 

^^^^ uit) = f{u{t),t), te{o,n 

w(0) = uo, 

where u : [0, T] is the solution to be computed, uq e a given initial 

condition, T > a given final time, and / : x (0, T] a given function that 

is Lipschitz-continuous in u and bounded. 

In the previous two parts of our series on multiadaptive Galerkin methods, we 
proved a posteriori error estimates, through which the time steps are adaptively de- 
termined from residual feedback and stability information, obtained by solving a dual 
linearized problem. In this paper, we prove a priori error estimates for mcG{q) and 
mdG((7). We also prove the stability estimates and interpolation estimates which are 
essential to the a priori error analysis. 

Standard methods for the time-discretization of ()1.1|) require that the resolution 
is equal for all components Ui{t) of the computed approximate solution U{t) of (jl.ip . 
This includes all standard Galerkin or Runge-Kutta methods; see |4l [23l [24l |4TJ 
[2]. Using the same time step sequence k — k{t) for all components could become 
very costly if the different components of the solution exhibit multiple time scales of 
different magnitudes. We therefore propose a new representation of the solution in 
which the difference in time scales is reflected in the componentwise time-discretization 
of (|l.ip . that is, each component Ui(t) is computed using an individual time step 
sequence ki = ki{t). 
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The multiadaptive Galerkin methods mcG{q) and mdG((7) first presented in [3 2) 
are formulated as extensions of the standard continuous and discontinuous Galerkin 
methods cG{q) and dG((7), studied earlier in detail by Hulme [551 HZ]; Jamet 
Delfour, Hager, and Trochu 7 , Eriksson, Johnson, and Thomee ^T6[ l30l [TTl [T2l [TOl 
[131 [H [m [g, and Estep et' al. [HI [H [IS [H [20]. As such, the analysis of the 
mcG{q) and mdG(g) methods can be carried out within the existing framework, but 
the extension to multiadaptive time-stepping leads to some technical challenges, in 
particular, proving the appropriate interpolation estimates. 

Local (multiadaptive) time-stepping has been explored before to some extent for 
specific applications, including specialized integrators for the n-body problem [37l[5l[T] 
and low-order methods for conservation laws [311 1211 [S]- Early attempts at local 
time-stepping include [IS] [25]. Recently, a new class of related methods, known as 
asynchronous variational integrators (AVI) with local time steps, has been proposed 

m- 

1.1. Main results. The main results of this paper are a priori error estimates 
for the mcG{q) and mdG((7) methods, respectively, of the form 

(1-2) ll^(r)lk.<C5(T)||fc^^u(2^)|L^([o,^],^) 

and 

(1-3) Mnw < ^^^ml|fc"'+'"'"'+''IL^([o,T],;,) 

for p = 2 or p — oo, where C is an interpolation constant, S{T) is a (computable) sta- 
bility factor, and /c^'^w^^'^ (or k'^i'^^u^^'^'^^'') combines local time steps ki = ki{t) with 
derivatives of the exact solution u. The norm L^{I, || • ||) is defined by ||w||i^(/n.||) = 
supjg/ ||w(i)||. These estimates state that the mcG(g) method is of order 2q and that 
the mdG{q) method is of order 2q + 1 in the local time step. We refer to section 
16.21 for the exact results. It should be noted that superconvergence is obtained only 
at synchronized time levels, such as the end-point t = T . For the general nonlinear 
problem, we obtain exponential estimates for the stability factor S{T). In [S^, we 
prove that for a parabolic model problem, the stability factor remains bounded and 
of unit size, independent of T (up to a logarithmic factor) . 

1.2. Notation. The following notation is used throughout this paper. Each 
component Ui{t), i — 1, . . . , N, the approximate m{c/d)G{q) solution U{t) of (II. ip 
is a piecewise polynomial on a partition of (0,T] into Mi subintervals. Subinterval j 
for component i is denoted by lij — {tij-i,tij], and the length of the subinterval is 
given by the local time step kij — tij — tij_i. This is illustrated in Figure[Tl On each 
subinterval lij, Ui\i-j is a polynomial of degree qij and we refer to (lij, Ui\i.^) as an 
element. 

Furthermore, we shall assume that the interval (0, T] is partitioned into blocks 
between certain synchronized time levels = To <?!<••• < Tm — T . We refer to 
the set of intervals Tn between two synchronized time levels r„_i and T„ as a time 
slab: 

7n — • ^71 — 1 ^ tij — i <^ tij ^ Tji^ . 

We denote the length of a time slab by Kn — Tn — Tn-i- We also refer to the entire 
collection of intervals lij as the partition T ■ 

Since different components use different time steps, a local interval lij may contain 
nodal points for other components, that is, some tiiji e (ti,j^i,tij). We denote the 
set of such internal nodes on a local interval by Afij . 
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Fig. 1. Individual partitions of the interval (0,T] for different components. Elements between 
common synchronized time levels are organized in time slabs. In this example, we have N = 6 and 
M = 4. 



1.3. Outline of the paper. The outline of this paper is as follows. In section[2l 
we give the full definition of the multiadaptive Galerkin methods mcG{q) and mdG{q). 
We also introduce the dual methods mcG{q)* and mdG((7)*, which are of importance 
to the a priori error analysis. In sections[3]and[4l respectively, we then prove existence 
and stability of the discrete solutions as defined in section [2l 

In section [5l we prove the interpolation estimates that we later use to prove the 
a priori error estimates in section [6l Proving the interpolation estimates is technically 
challenging, since the function to be interpolated may be discontinuous within the 
interval of interpolation. To measure the regularity of the interpolated function, it is 
then necessary to take into consideration the size of the jump in function value and 
derivatives at each point of discontinuity. 

Finally, in section [71 we present some numerical evidence for the a priori error 
estimates by solving a simple model problem and showing that we obtain the predicted 
convergence rates, k'^'^ and respectively, for the mcG{q) and mdG((7) methods. 

2. Definition of methods. In this section, we give the definitions of the mul- 
tiadaptive Galerkin methods mcG{q) and mdG((7). The multiadaptive methods are 
obtained as extensions of the standard (monoadaptive) Galerkin methods cG{q) and 
dG(g) by extending the trial and test spaces to allow individual time step sequences 
for different components. 

As an important tool for the a priori error analysis in section [6l we also introduce 
the discrete dual problem and the discrete dual methods mcG{q)* and mdG((7)*. 

2.1. Multiadaptive continuous Galerkin, mcG(q). To formulate the mcG((7) 
method, we define the trial .space V and the test .space V as 

y = {«e [C([0,T])]^ :^;,|,.^ gT''- (/„■), J = 1, . . . , M„ i = 1, . . . , TV}, 
V={v:v,\u^ ^V^^'-^h,), j = l,...,M„ z=l,...,7V}, 

where T"^{I) denotes the linear space of polynomials of degree q on an interval / C M. 
In other words, V is the space of vector- valued continuous piecewisc polynomials of 
degree q = {qi{t)) with qi{t) > 1 on the partition T, and V is the space of vector- 
valued (possibly discontinuous) piecewise polynomials of degree q — 1 = {qi(t) — 1) on 
the same partition. 
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We now define the nicG(g) method for (jl.ip as follows: Find U with U{0) = 
uq such that 

(2.2) / {U,v)dt= [ {f{U,-),v)dt \fveV, 
Jo Jo 

where (•, •) denotes the inner product. With a suitable choice of test function v, it 
follows that the global problem (|2.2[) can be restated as a sequence of successive local 
problems for each component: For i — 1, . . . , N , j = 1, . . . , Mi, find Ui\i.- G V^^' {hj) 
with Ui{ti_j^i) given such that 

(2.3) / U^vdt= I h{U,-)vdt yveV'^^-^ihj), 



where the initial condition is specified for j = 1, . . . , by Ui{0) — Ui{0). 

We define the residual R of the approximate solution U by Ri{U,t) — Ui{t) 
fi{U{t),t). In terms of the residual, we can rewrite (12. 3p in the form 

(2.4) / R,iU,-)vdt^O yv eV'^^-^ihj), j = l,...,M„ i^l,...,N, 



that is, the residual is orthogonal to the test space on each local interval. We refer to 
(|2.4p as the Galerkin orthogonality of the mcG((7) method. 

2.2. Multiadaptive discontinuous Galerkin, mdG(q). For mdG(g), we de- 
fine the trial and test spaces by 

(2.5) V = {v:v,\u^€V'^^'{h,), j = 1, . . . , M„ z = 1, . . . , 7V}, 

that is, both trial and test functions are vector-valued (possibly discontinuous) piece- 
wise polynomials of degree q = {qi{t)) with qi{t) > on the partition T. By definition, 
the mdG(g') solution U Q V is left-continuous. 

We now define the mdG{q) method for (jl.ip as follows: Find U ^ V with U{0^) = 
Uq such that 



N Mi 

(2-6) EE 

i=l j=l 



/ iliVi dt 


-I 







i-T 

{f{U,-),v)dt VveV, 



where [Ui]ij-i — Ui{t'l_j_^) — Ui{i^^j-i) denotes the jump in Ui{t) across the node 
t = iij-i, and where v(t~'') = lim5_^(+ v{s). 

The mdG(g) method in local form, corresponding to (|2.3p . reads as follows: For 
i = 1, . . . , A^, j = 1, . . . , M,, find f/,|/,^. G V'^ {h-j) such that 

(2.7) [C/,],,,_i«(t,.j_i) / U,vdt= [ MU,-)vdt Vz;eP«- (/,,), 

where the initial condition is specified for i = 1, . . . , iV by Ui{Q~) — Ui{Q). 

The residual R is defined on the inner of each local interval lij by Ri(U,t) = 
Ui{t) — fi{U{t),t). In terms of the residual, p.7p can be restated in the form 

(2.8) [U,],,,.iv{t+^_^)+ f R,iU,-)vdt^O yveV'^'ihj) 

J iij 

for j = 1, . . . , Mi, i = 1, . . . , N . We refer to (12. 8p as the Galerkin orthogonality of the 
mdG((7) method. 
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2.3. The dual problem. The dual problem is the standard tool for error anal- 
ysis, a priori or a posteriori, of Galerkin finite element methods for the numerical 
solution of differential equations; see [Sl[3]. For the a posteriori error analysis of the 
niultiadaptive Galerkin methods mcG(g) and TadG{q) in [32], we formulate a continu- 
ous dual problem. For the a priori error analysis of this paper, we formulate instead a 
discrete dual problem. The discrete dual problem was first introduced for the family 
of discontinuous Galerkin methods dG((7) in As we shall see, the discrete dual 
problem can be expressed as a Galerkin method for a continuous problem. 

The discrete dual solution $ : [0,T] — >■ M.^ is a Galerkin approximation of the 
exact solution : [0,r] — ?> M.^ of the continuous dual problem 

~<Pit)^J^{7ru,U,t)(b{t)+g{t), te[0,T), 
0(r) = ^, 

where nu is an interpolant or a projection of the exact solution u of (II. ip . g : [0, T] — > 
is a given function, ip G is a given initial condition, and 

(2.10) J^{Tru,U,t) ^ ( [ ^{snu{t) + {l- s)U{t),t)ds 



du 

that is, an appropriate mean value of the transpose of the Jacobian of the right-hand 
side f{-,t) evaluated at TTu{t) and U{t). Note that by the chain rule, we have 

(2.11) Jinu, U, ■)iU - nu) = /(C/, •) - /(ttu, •)■ 

The data (■0, g) of the dual problem allow us to obtain error estimates for different 
functionals L^^g of the error e — U — u. 

We define below two new Galerkin methods for the dual problem (|2.9p : the dual 
methods mcG(g)* and nidG(g)*. We will later use the mcG{q)* method to express 
the error of the mcG{q) solution of (jl.ip in terms of the mcG{q)* solution of (|2.9p . 
Similarly, we will express the error of the mdG(g) solution of (II. ip in terms of the 
nidG(g)* solution of (f^ . 

2.4. Multiadaptive dual continuous Galerkin, mcG(q)*. In the formula- 
tion of the dual method of mcG{q), we interchange the trial and test spaces of nicG((7). 
With the same definitions of V and V as in (|2.ip . we thus define the mcG(g)* method 
for dMl as follows: Find $ G V with $(T+) = V such that 

(2.12) / {v,^)dt^ [ {J{nu,U,-)v,^)+L^Jv) 
Jo Jo 

for all u e y with v{0) — 0, where 

(2.13) L^.g{v) = {v{T),ij)+ [ {v,g)dt. 

Jo 

Notice the extra condition that the test functions should vanish at t = 0, which is 
introduced to make the dimension of the test space equal to the dimension of the trial 
space. Integrating by parts, p.l2p can alternatively be expressed in the form 



N Ah 

(2-14) EE 

i=l 3=1 L 



Jlii Jo 



i-T 

{J^{7Tu,U,-)<i> + g,v) dt. 
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2.5. Multiadaptive dual discontinuous Galerkin, mdG(q)*. With the same 
definitions of V and V as in (I2.5p . we define the mdG(g)* method for (|2.9p as follows: 
Find $ e with $(r+) = V such that 



N Mi 

(2.15) EE 

2 = 1 j=l 



/ 'i'l^'i dt 


-I 







(J(7ru, f/, ■)v, $) dt + i^,g(w) 



for aWveV with t;(0 ) = 0. Integrating by parts, (|2.15p can alternatively be expressed 
in the form 



(2.16) 



N Mi 

EE 

i=i ]=i 



<^,v, dt 



[J^ {■Ku,U,-)<^ + g,v)dt. 



3. Existence of solutions. To prove existence of the discrete mcG(g), mdG((7), 
mcG(<7)*, and TcidG{q)* solutions defined in the previous section, we formulate fixed 
point iterations for the construction of solutions. Existence then follows from the 
Banach fixed point theorem if the time steps are sufficiently small. 

Lemma 3.1 (fixed point iteration). Let Tn be a time slab with synchronized 
time levels T„_i and Tn. With time reversed for the dual methods (to simplify the 
notation), the m.cG{q), mdG(g), mcG(g)*, anduidG{q)* methods can all be expressed 
in the following form: For all lij G Tn, find {^ijn} (the degrees of freedom for Ui on 
lij ) such that 



(3.1) 



^ij7l = Ui{0) + 



ti, j-i 



fr{U,-)dt+ / w[f-l(n,(<))/,(C/,.)di, 



where Tij{t) = (t — tij^i)/(tij — tij^i) and {wn'^'} is a set of polynomial weight 
functions on [0, 1] . 

Proof. The result follows from the definitions of the mcG{q), mdG{q), nicG{q)* , 
and mdG((7)* methods, using an appropriate basis for the trial and test spaces. See 
P] for details. □ 

Theorem 3.2 (existence of solutions). Let I'C = max i4r„ be the maximum time 
slab length and define the Lipschitz constant Lf > by 



(3.2) 
// now 
(3.3) 



\\f{x,t)-f{y,t)\\i^<Lf\\x-y\\i^ Vie [0,r] Vx, y G 



KCLf < 1 



where C — C{q) > is a constant depending only on the order and method, the 
fixed point iteration p.l|) converges to the unique solution of (|2.2|) . (|2.6|) . ()2.12|) . and 
(j2.15|) . respectively. 

Proof. The result follows by Lemma 13.11 and an application of the Banach fixed 
point theorem. See [Mi for details. □ 

4. Stability of solutions. Write the dual problem (|2.9p for cj) — (f){t) in the 



form 
(4.1) 



0(T) = ^. 
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For simplicity, we consider only the case g = 0. With = (p{T — t), we have 

■w{t) = —(j){T — t) ^ —A^{T — t)w{t), and so (|4.1I) can be written as a forward 
problem for w in the form 

wit) + B{t)w{t)^o, te{o,n 

'w(0) = Wo, 

where wq = t/j and = {T — t). Below, w represents either u or (/)(T — •) and, 
correspondingly, W represents either the discrete mc/dG(g) approximation U of u or 
the discrete mc/dG(<7)* approximation $ of (p. 

4.1. A general exponential estimate. The general exponential stability esti- 
mate is based on the following version of the discrete Gronwall inequality. 

Lemma 4.1 (discrete Gronwall inequality). Assume that z,a : N — > M are non- 
negative, a(m) < 1/2 for all m, and z{n) < C + '^ln^ia{m)z{m) for all n. Then 
z{n) < 2Cexp(E::,^\ 2a(m)) for n = 1,2, ... . 

Proof. By a standard discrete Gronwall inequality [3S], z{n) < C exp{J2l^^J^^ a{m)) 
if z{n) < C + X]m=o '^("^)-^('^) for n > 1 and z(0) < C. Here, (1 — a{n))z{n) < 
C + X]m=i a{m)z{m), and so z{n) < 2C + X]m=i '^0'{m)z{m), since 1 — a{n) > 1/2. 
The result now follows if we take a(0) = 2:(0) = 0. □ 

Theorem 4.2 (stability estimate). Let W be the mcG(9), mdG(g), mcG(g)*, or 
mdG((3')* solution of (j4.2l) . Then there is a constant C — C{q), depending only on 
the highest order maxg^, such that */ ^nC'||5||L^([T„_i,T„],;p) ^ 1 for n — 1, . . . ,M, 
then 

(4-3) ||W^||l^([T,._i,T„]J,) < C\\wq\\i^ exp [ ^ ^mC^I|-B||L„([T„_i,T,„],i,) J 

for n — 1, . . . , M , 1 < p < oo. 

Proof. By Lemma l3Tl we can write the mcG{q), mdG(g), mcG(q)*, and nidG((7)* 
methods in the form ^ w,{0) + J^;'-'-' f^iW, ■) dt + J^^^ w^^r\n,{t))f,{W, ■) dt. 
Applied to the linear model problem (|4.2I) . we have S^ijn' — Wj(0) — /q''^~^ {BW)i dt— 
!i.w^n''\n3{i)){BW)^ dt, and so 



ICijn'l < \Wi 



{BW)^ dt 



w^:^l^\n,{t)){BW\dt 



< \w^{0)\ +C \iBW),\dt < \w^iO)\ + C /" " \iBW),\dt, 
Jo JO 



for aU t e [T„_i,T„], we have \W^{t)\ < C\w,{0)\ + C J^^ \{BW).^\ dt, and so 
\\W{t)\\i^ < C\\wo\\i^ + C \\BW\\i^ dt = C\\wo\\i^ +CJ2 II^W^II^P dt. 

Jo 1 Jt^-1 



The result now follows by letting Wn = ||T^||L^([T„_i,T„],ip)- D 

Remark 4.1. See |34| for an extension to multiadaptive time-stepping of the 
strong stability estimate Lemma 6.1 for parabolic problems in [11) . 
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5. Interpolation estimates. In this section, we introduce a pair of carefully 
chosen interpolants, 7r|.Q and tt^Jq, which are central to the a priori error analysis 



of the mcG{q) and mdG(g) methods. The interpolants are defined in section 15.11 
In section 15.21 we discuss the interpolation of piecewise smooth functions, that is, 
the interpolation of functions which may be discontinuous within the interval of in- 
terpolation, and then present the basic general interpolation estimates for the two 
interpolants 7r|,Q and tt^q. 

For the a priori error analysis of the mcG{q) and mdG((7) methods, we will also 
need a special interpolation estimate for the function if — J^$, where J is the 
Jacobian of the right-hand side / of (jl.ll) and $ is the discrete dual solution as 
defined in section [51 including estimates for the size of the jump in function value and 
derivatives for the function ip at points of discontinuity. These estimates are proved in 
section [5751 based on a representation formula for the mcG(g) and mdG(g) solutions 
of ([TT]) . 



5.1. Interpolants. The interpolant 7r|,Q : V 
lowing conditions: 



P9([a,6]) is defined by the fol- 



(5.1) 



7rt«lw(a) = v{a) and 7r[|i;(fe) = v{b), 



where V denotes the set of functions that are piecewise and bounded on [a,b]. 
In other words, 7r|,QW is the polynomial of degree q that interpolates v at the end- 
points of the interval [a, b] and additionally satisfies q — 1 projection conditions. This 
is illustrated in Figure [51 We also define the dual interpolant 7r|^Q, as the standard 
L2-projcction onto T"^~^{[a, b]). 




0.1 0.2 0.3 



.2 0.3 0.4 



0.6 0.7 O.a 9 



Fig. 2. The interpolant tt^q"^ (dashed) of the function v(x) = x sin(7a;) (solid) on [0,1] for 
g = 1 (left) and q = 3 (right). 

The interpolant tt^q : 1/ — > 7^^([a, 6]) is defined by the fohowing conditions: 



(5.2) 



'dG 

4t«(6)=T;(5), 



Qomial of degree q that interpolates v at tl 
the interval [a, b] and additionally satisfies q projection conditions. This is illustrated 



that is, t^I^qV is the polynomial of degree q that interpolates v at the right end-point of 
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Fig. 3. The interpolant t^^qV (dashed) of the function v{x) = x sin(7x) (solid) on [0,1] for 
(7 = (left) and q = 3 (right). 




Fig. 4. A piecewise smooth function v and its interpolant nv. 



in Figure [3l The dual interpolant tt^jq, is defined similarly, with the difference being 
that the left end-point x — a is used for interpolation. 

5.2. Basic interpolation estimates. To estimate the size of the interpolation 
error ttv — v for a given function v, we express the interpolation error in terms of the 
regularity of v and the length of the interpolation interval, k = b — a. Specifically, 
when V E ^^^([a, b]) C V for some g > 0, we obtain estimates of the form 

(5.3) ||(7rw)(f) < C/fc«+i-f||z;(«+i)||, p = 0, . . . , g + 1, 

where |1 • |1 = || • |lL^([a.fc]) denotes the maximum norm on [a,b]. This estimate is 
a simple consequence of the Peano kernel theorem [40^ if one can show that the 
interpolant tt : — )■ ^^([a, 6]) C V is linear and bounded on V and that tt is exact on 
Viia^b]) C V, that is, TTu = u for all v G 7'«([a,6]). 

In the general case, where the interpolated function v is only piecewise smooth 
(see Figure HI, we also need to include the size of the jump in function value 

and derivatives at each point x of discontinuity within (a, b) to measure the regularity 
of the interpolated function v. In [34j, we prove the following extensions of the basic 
estimate (|5.3|) . 

Lemma 5.1. If n is linear and bounded on V and is exact on T"^{[a,b]) C V, 
then there is a constant C — C{q) > such that for all v piecewise C'^^ on [a, b] with 
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Fig. 5. // some other component I ^ i has a node within lij, then <i>; may be discontinuous 
within lij, causing ipi to be discontinuous within lij. 

discontinuities at a < Xi < ■ ■ ■ < Xn < b, 

n r 

(5.4) ||(7ru)(P) - < CF+i-^'||w(''+i)|| ^ ^ | 

j = l m=0 

for p — 0, . . . ,r + 1, r = 0, . . . , q. 

Lemma 5.2. If n is linear and bounded on V and is exact on 7"([a, 6]) C V, 
then there is a constant C — C(q) > such that for all v piecewise C"*"^ on [a, b] with 
discontinuities at a < xi < ■ ■ ■ < Xn < b, 

n p—l 

(5.5) IKttw)^"'!! < cy^^w + cY^ 

j — 1 ni—O 

for p ^ 0, . . . ,q. 

Lemmas 15.11 and 15.21 apply to both the 7r|.Q interpolant (for q > 1) and the tt^q 
interpolant (for q > 0) defined in sectfon lOl The hnearity of both interpolants fohows 
directly from the definition of the interpolants. The proofs that both interpolants are 
bounded and exact on V^a, b]) are given in detail in [34 and [35 . 

5.3. A special interpolation estimate. To prove a priori error estimates for 
mcG(g) and mdG(g) in section |6l we need to estimate the interpolation error mp — ip 
for the function ip defined by 

AT 

(5.6) (p^ = (J^(7rM,u, •)*)»= XI « = 

1=1 

We note that Lpi may be discontinuous within if contains a node for some other 
component, which is generally the case. This is illustrated in Figure [SJ Note that on 
the right-hand side / is linearized around a mean value of -ku and u. 

An interpolation estimate for inp—Lp follows directly from Lemma l5.1l To use this 
estimate, we need to estimate the size of the jump in function value and derivatives at 
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each internal node tij of the partition T. To obtain this estimate, we need to make a 
number of additional assumptions on the right-hand side / of (|l.ip and the partition 
T. These assumptions are discussed in section 15.3.21 Based on the assumptions and 
the representation formula presented in section 15.3.11 we obtain the jump estimates 
in section [5.3.31 and, finally, in section [5.3.41 the interpolation estimate for ip. 

5.3.1. A representation formula. The proof of jump estimates for the multi- 
adaptive Galerkin methods mcG{q) and mdG(g) is based on expressing the solutions 
as certain interpolants. These representations are obtained as follows. Let U be the 
mcG{q) or mdG{q) solution of (jl.ll) and define, for i = 1, . . . , N , 

(5.7) U,{t)^u,{0)+ f MUis),s)ds. 

Jo 

Similarly, for $ the mcG{q)* or mdG{q)* solution of (|2.9p . we define, for i = 1, . . . , A'^, 

(5.8) $,(*)= V^,+^ f*ms),s)ds, 

where /*($, •) = {nu, U, •)$ + g. We note that U = /{U, •) and -$ = /*($, •). 

It now turns out that U can be expressed as an interpolant of U. Similarly, $ 
can be expressed as an interpolant of <&. We present these representations in Lemmas 
15.31 and 15.41 We remind the reader about the interpolants 7r|,Q, 7r|.Q,, tt^J^q, and tt^j^q. , 
defined in section [STTl 

Lemma 5.3. The mcG(g) solution U of (11. ip can expressed in the form U = 7r|.Qt7. 
Similarly, the mcG((7)* solution $ of (|2.9p can be expressed in the form $ — 7r|^Q,$, 

that is, Ui = Trj.Q^'c/i and $i = 7r|,Qi'$i on each local interval Iij. 

Proof. The representation formulas follow by the definitions of the mcG(g) and 
mcG(g)* methods and the interpolants 7r|,Q and 7r|,Q,. See [34 for details. □ 

Lemma 5.4. The mdG(q) solution U of (jl.ll) can expressed in the form U — t^I^qU- 
Similarly, the mdG(g)* solution $ of (|2.9p can be expressed in the form $ — ttI^q*^, 

that is, Ui = TT^jQ^'cA and $i = T^^^Q}'^i on each local interval lij. 

Proof. The representation formulas follow by the definitions of the mdG(q) and 
mdG((3')* methods and the interpolants tt^j^q and tt^^q, . See [34] for details. □ 

5.3.2. Assumptions. To estimate the size of the jump in function value and 
derivatives for the function ip defined in (j5.6p . we make the following assumptions. 
Given a time slab T, assume that for each pair of local intervals lij and within 
the time slab, we have 

(Al) qij = q„in = q 

and 

(A2) kij > a kmn 

for some q > and some a G (0, 1). The dependence on a in the error estimates is 
weak (see Remark 15. 11) . so assumption (A2) does not prevent multiadaptivity. 
We also assume that the problem (jl.ip is autonomous, 

(A3) dfJdt^O, i^l,...,N, 
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noting that the dual problem nevertheless will be nonautonomous in general. Fur- 
thermore, we assume that 

(A4) ||/*||c«+i(r) < oo, i^l,...,N, 

where || • ||_dp(7-) is defined for v : M.^ — > M and p > by ||f||£ip(T) = rnax„=o,...,p 
ll-^"^IUoo(T,;oo)' '^it^ the norm \\D'^v\\L^(^r.i^) defined by IjZ?"^ ••• 1,^(7-) < 
\\D''v\\lZ{t[i2)\\w'^\\i^ for a£ . . , w"" G K^, and the nth-order 

tensor given by 



dxi^ ■ ■ ■ dx 



ii — l in—1 

Furthermore, we choose C/ > maxi=i^..._Ar ||/i||_Dfl+i(r) such that 

(5.9) \\dP/dtP{df/duY{xm\i^ <CfCP 

for p = 0, . . . , q, and 



(5.10) \\[dP/dtP{df/duyixitml^ < Qf]cr"||[a;(") 



tlUo 



n=0 

pN 



for p = 0, . . . , g — 1 and any given a; : M — > R , where Cx > denotes a constant such 
that ll^;*-"' ||l^(7-,;^) < C" for n = 1, . . . ,p. Note that C/ = Cf{t) defines a piecewise 
constant function on the partition = Tq < Ti < • • • < Tm — T. Note also that 
assumption (A4) implies that each is bounded by C/. 

We further assume that there is a constant > such that 

(A5) kijCf < Ck 

for each local interval lij. We summarize the list of assumptions as follows: 
(Al) the local orders qij are equal within each time slab; 
(A2) the local time steps kij are semiuniform within each time slab; 
(A3) / is autonomous; 
(A4) / and its derivatives are bounded; 
(A5) the local time steps kij are small. 

5.3.3. Estimates of derivatives and jumps. To estimate higher-order deriva- 
tives, we face the problem of taking higher-order derivatives of f{U{t), t) with respect 
to t. In Lemmas 15.51 and 15.61 we present basic estimates for composite functions vox 
with V : — R and a; : R — > R^. The proofs are based on a straightforward 
application of the chain rule and Leibniz rule and are given in full detail in ,34. . 

Lemma 5.5. Let v : R-'^ — > R 6e p > times differentiable in all its variables, 
let X : M. R^ be p times differentiable, and let Cx > be a constant such that 
ll^'""''IUoo(K,icxj) — f'^''" = 1; ■ • ■ )P- Then there is a constant C = C{p) > such 
that 

(5.11) 



dP{v o x) 



dtp 



Lemma 5.6. Let v : M.^ — >R&ep+l>l times differentiable in all its variables, 
Zei X : R — > R^ be p times differentiable, except possibly at some t G R, and let 
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Cx > be a constant such that ||a;''"-'||L<^(R,io„) < C^x f^^ ^ ~ 1; • ■ • iP- Then there is 
a constant C = C(ja) > such that 

<cii-iiD.+H«)E^rllMJL^. 

n=0 

We now prove estimates for derivatives and jumps of the mcG(g) or mdG(g) 
solution U of the general nonlinear problem under the assumptions listed in 

section [5.3.2l Similarly, one can obtain estimates for the discrete dual solution $ and 
the function ip defined in (j5.6p . from which the desired interpolation estimates follow. 

To obtain estimates for the multiadaptive solution U, we first prove estimates for 
the function U defined in section [5.3. II The estimates for U then follow by induction. 

To simplify the estimates, we introduce the following notation. For given p > 0, 
let Cu,p > C/ be a constant such that 

(5.13) ||f/(")||^^(^,^) <C^_p, n = l,...,p. 

For p = 0, we define Cu,o = Cf. Temporarily, we assume that there is a constant 
c'j, > such that for each p, 

(A5') h,Cu,p < 4. 

This assumption will be removed in Lemma 15.91 In the following lemma, we use 
assumptions (Al), (A3), and (A4) to derive estimates for U in terms of Cu,p and C/. 

Lemma 5.7 (derivative and jump estimates for U). Let U be the mcG(q) or 
mdG{q) solution of and define U as in (|5.7|) . If assumptions (Al), (A3), and 

(A4) hold, then there is a constant C — C{q) > such that 

(5.14) ||C/(P)||^^(^^^^)<CC^_P_1, p = l,...,9-+l, 
and 

(5-15) l|[^^^^]*_JL„<^^E^rp-ill[C^^"^]*_JL.' p=l,...,g-+l, 

ra=0 

for each local interval lij , where iij-i is an internal node of the time slab T. 

Proof. By definition, U^^'^ = ^^p-i fi{U), and so the results follow directly by 
Lemmas 15.51 and 15.61 noting that C/ < Cu,p^i- □ 

By Lemma 15.71 we now obtain the following estimate for the size of the jump in 
function value and derivatives for U. 

Lemma 5.8 (jump estimates for U). Let U be the mcG(g) or mdG{q) solution 
of (jl.ip . // assumptions (A1)-(A5) and (A5') hold, then there is a constant C = 
C{q,Ck,c'f.,a) > such that 

(5.16) \\[U^''^u,.J\i^<Cki;'-''Cl+\ p = 0,...,r + l, r = 0,...,g, 

for each local interval lij , where iij-i is an internal node of the time slab T- 
Proof. The proof is by induction. We first note that at t — tij-i, we have 

[u^\ = u^'Ht^) - ul'\t^) + u^'\t+) u^^\t-) + ul^\t-) ul^\t~) 

= e+ + eo + e_. 



(5.12) 



d^{v o x) 
dtp 
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By Lemma 15.31 (or Lemma 15 .4^ . U is an interpolant of U and so, by Lemma |5.1[ we 
have 

xeN'ij m=l 

for p — 0, . . . ,r and r — 0, . . . , q. Note that the second sum starts at m = 1 rather 
than at m = 0, since U is continuous. Similarly, we have 

x£N'i, j-1 m=l 

To estimate eo, we note that cq — for p — 0, since U is continuous. For p = 1, . . . , q+ 

l,LemmaOgives|eo| = \[U^%\ < C Ef.^o C&7p" i II [f/("^]t Ihoo • By assumption (A2), 
it then follows that (|5.16p holds for r = 0. 

Assume now that (|5.16p holds for r = f — 1 > 0. Then, by Lemma 15.71 and 
assumption (A5'), it follows that 

f m — 1 

ie+i < cfc^+i-'^c^+i + ^ E E >'T' E Q™" ilim^L 

x^Afij m—l n— 

<r nh^+'^-Pnf+l I im-pf-mi-n , (r'-l) + l-n^(f-l) + l 



U,f ■ 



Similarly, we obtain the estimate |e_| < Ck^^^ ^C\j^^. Finally, we use Lemma [5.71 
and assumption (A5') to obtain the estimate 

leoi < cY^cfj:;_,\\[u-],i^ < cj2c!i:;iikt'^^'-''c^^T 

p-1 
n=0 

Summing up, we thus obtain \[U^%\ < \e+\ + |eo| + |e_| < Ckl+^~^Cl+\ and so 
(|5.16p follows by induction. □ 

By Lemmas 15.71 and 15.81 '^^ now obtain the following estimate for derivatives of 
the solution U. 

Lemma 5.9 (derivative estimates for t/). LetU be themcG{q) orindG{q) solution 
of (jl.ip . If assumptions (A1)-(A5) hold, then there is a constant C = C{q,Ck,a) > 
such that 

(5.17) ||f/(^)||,^(^_,^^<CC,^ p^l....,q. 

Proof. By Lemma l5.3l for Lemma l5.4p . U is an interpolant of U and so, by Lemma 
15.11 we have 

x^Afij m—l 
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for some constant C = C'{q). For p — 1, we thus obtain the estimate 

mh^iU,) < C'||L^.|U^(/,,) = C'||/.(C/)|U^(,.^) < C'Cf 

by assumption (A4), and so (|5.17l) holds for p = 1. 

For p = 2, . . . ,q, assuming that (A5') holds for Cu,p-i, we use Lemmas 15.71 and 
(with r = p — 1) and assumption (A2) to obtain 



p— 1 m—l 

where C = C{q, Ck, c'l^, a). This holds for all components i and all local intervals lij 
within the time slab T, and so 

where by definition Cu,p~i is a constant such that ||J/*-"''||L^(T,i3o) — ^Up~i ^'^^ = 
1, ... ,p — 1. Starting at p = 1, we now define Cu,i = C'lC'/ with Ci = C" = C'{q). It 
then follows that (A5') holds for C(7,i with c'j, = C'c^, and thus 

ll^*'^^llLoo(r,ioo) - ^^u.2-i = C'C'&a = (^2^1, 

where C2 = C2{q,Ck,a). We may thus define Cu,2 = max(CiC/, \fC2Cf ). Continu- 
ing, we note that (A5') holds for Cu,2^ and thus 

11^^'^^ llioo(T,ioo) — ^^u,?,-i — CC^2 = C'sC'/, 
where C3 — C^{q^ c^, a). In this way, we obtain a sequence of constants Ci, . . . , Cg, 
depending only on g, c^, and a, such that ||C^'''-'||L^(T,i3o) — ^p^^f for p = 1, . . . , 
and so (|5.17l) follows if we take C = max^^i^. .^^ Ci. □ 

Having now removed the additional assumption (A5'), we obtain the following 
version of Lemma 15.81 

Lemma 5.10 (jump estimates for U). Let U he the mcG(g) or mdG(g) solution 
of (jl.ip . If assumptions (Al)-(A5) hold, then there is a constant C = C(g, Cfc,a) > 
such that 

(5.18) \\[U'^'^]u,,JL<Ck!;'-^Cj+\ p = 0,...,g- 

for each local interval lij, where ijj-i is an internal node of the time slab T. 

Similarly, we obtain estimates for the discrete dual solution $ and the function 
if. In Lemma 15.111 we present the estimates for the function ip. 

Lemma 5.11 (estimates for ip). Let ip be defined as in (j5.6p . // assumptions 
(A1)-(A5) hold, then there is a constant C =^ C{q,Ck,a) > such that 

(5-19) lk!'^IL^(/.,)<^^^r'll*llio^(r,^o.), P = 0,...,g.„ 

and 

(5.20) \[^'f^]J<Ck:;^-''cy'^'mL^(^rA^) VxeM,, p = o,...,<z.,-i, 

with rij = qij for the nicG{q) method and rij = qij + 1 for the mdG(g) method. This 
holds for each local interval Lij within the time slab T. 
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5.3.4. Interpolation estimates. Using the basic interpolation estimate of sec- 
tion 15.21 we now obtain the following important interpolation estimates for the func- 
tion Lp. 

Lemma 5.12 (interpolation estimates for Lp). Let Lp be defined as in (|5.6I) . // 
assumptions (Al)-(A5) hold, then there is a constant C — C{q,Ck,a) > such that 

(5.21) ||7r[^(^^-'V.-¥'»|L^(,^^) <Cfc^''''Cf ll<i>llL^(r.;„.), q^J^q>2, 
and 

(5.22) ||4g~'V. - ¥'.|L^(,,^.) < CkT/Cf^^'mL^ir.i^), q^j =q>l, 

for each local interval lij within the time slab T- 

Proof. To prove (j5.2ip . we use Lemma TS-li with r = qij — 2 and p — 0, together 
with Lemma 15.111 to obtain 

x£j\fij m=0 

x^Nij m=0 

from which the estimate follows. The estimate for -k^^^ ^^(pi — ipi is obtained simi- 
larly. □ 

Remark 5.1. Note that there is only a weak dependence on Ck and a, since the 
jump term contains an extra factor kij . If higher- order terms can be ignored, then the 
dependence on Ck and a can be removed. 

6. A priori error estimates. To prove a priori error estimates for the v[\cG{q) 
and mdG(g) methods, we derive error representations in section [6.11 and then obtain 
the a priori error estimates in section [6.21 for the general nonlinear case. We refer to 
|34| for a sharp a priori error estimate in the case of a parabolic model problem. 



6.1. Error representation. For each of the two methods, mcG(q) and mdG(g), 

we represent the error in terms of the discrete dual solution $ and an interpolant ttu of 

[<?] [9] 
t^qM or vru = ii^^u 



the exact solution u of (II. ip . using the special interpolants ttu = '"'cg''^ o^' ""'^ = 



defined in section [5j 

We write the error e = U — u \n the form 

(6.1) e = e -|- (ttu — u), 

where e = U — nu is represented in terms of the discrete dual solution and the residual 
of the interpolant. An estimate for the second part of the error, ttu — u, follows directly 
from an interpolation estimate. 

In Lemma 16. 1[ we present the error representation for the mcG(q) method, and 
then present the corresponding representation for the mdG(g) method in Lemma 16.21 
The error representations are obtained directly by choosing e as a test function for 
the discrete dual problems ()2.12p and ()2.15p . 

Lemma 6.1 (error representation for mcG(q)). Let U be the mcG(g) solution of 
(|1.1[) . let $ be the corresponding mcG(g)* solution of the dual problem (j2.9l) . and let 
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nu be any trial space approximation of the exact solution u of (|l.ip that interpolates 
u at the end-points of every local interval. Then 

L^^g{e) = {e{T),i,)+ j {e,g)dt^~( {R{^u, ■),^) dt, 
Jq Jo 

where e = U — ttu. 

Lemma 6.2 (error representation for mdG(g)). Let U be the nidG(g) solution of 
(jl.ll) . let $ be the corresponding mdG(g)* solution of the dual problem (|2.9p . and let 
t:u be any trial space approximation of the exact solution u of (|l.ip that interpolates 
u at the right end-point of every local interval. Then 

N Mi 

Li,,g{e) = -^Y^ 
i=i j=i 

where e = U — nu. 

With a special choice of interpolant, nu ~ t^]?qU and nu = '''IjqU, respectively, we 
obtain the following versions of the error representations. 

Corollary 6.3 (error representation for mcG(q)). Let U be the mcG(q) solution 
of (|l.ip and let $ be the corresponding mcG{q)* solution of the dual problem (j2.9p . 
Then 

Li>M^ f {f{niiu,-)-f{u,-),<S>)dt. 
Jo 

Proof. Integrate by parts and use the definition of the interpolant 7r|.Q. □ 
Corollary 6.4 (error representation for mdG(g)). Let U be the mdG(g) solution 

of (jl.ip and let $ be the corresponding mdG(g)* solution of the dual problem (j2.9p . 

Then 

L^.A^) = r •) - dt. 

Jo 

Proof. Integrate by parts and use the definition of the interpolant 7r|j^Q. □ 

6.2. A priori error estimates for the general nonlinear problem. Using 
the error representations of section [Ol the stability estimates of section |4l and the 
interpolation estimates of section [5l we now prove our main results: a priori error 
estimates for general order uYcG{q) and mdG((7). 

Theorem 6.5 (a priori error estimate for mcG((7)). Let U be the mcG(g) solution 
of (jl.ip and let $ be the corresponding mcG(g)* solution of the dual problem 
Then there is a constant C — C{q) > such that 



(6.2) \L^,gie)\<CS{T)\\k''+'u^'^+'^\ 



ioo([o,T],;2)' 



where {k'i+^u^i+^^)i{t) = k1'/^^\\u^^''^^''\\L^(u^) fort G 1,^ , and where the stability 

factor S(T) is given by S{T) = \\J^ (tt^^qU^u, dt. Furthermore, if assump- 

tions (Al)-(A5) hold, then there is a constant C = C{q, Ck, a) > such that 



(6.3) |Lv.,ff(e)| <C5(T)||fc295(29)| 



ioo([O.T],/i)' 
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where {k'^iu'^^i^)i{t) = k^^" Cf' ^||"i''''^^'||L„(/,,) for t G l,j , and where the stability 
factor S{T) is given by 

T M 

•^0 n=l 

Proof. By Corollary 16. 3[ we obtain 

L^,,(e) = r {f{nliu, •) - f{u, •),<!') dt = T (ttI^jU - u, {nl%u, u, •)$) d<. 
Jo Jo 

By Lemma [5. 1[ it now follows that 

<q|fc^+'s^+'||L^([o,T],/.) r||J^(7r[tu,w,.)$||,,dt, 



which proves (|6.2I) . To prove (|6.3p . we note that by definition, n^^Q^Ui—Ui is orthogonal 
to V'^'^^'^{Iij) for each local interval lij, and so, recalling that ip — {tt^^qU, u, •)$, 



where we take tt^^q "^^ipi = for qij — 1. By Lemmas 15.11 and 15.121 it now follows 

that 



where S{T) = /J' C/||$|U^(r,/^) = E^lii^f^nQll^llL^cr^.;^)- □ 

Similarly, we obtain the following a priori error estimate for the mdG(g) method. 
Theorem 6.6 (a priori error estimate for mdG(g)). Let U be the mdG(g) solution 

of (jl.ip and let $ be the corresponding mdG(g)* solution of the dual problem (|2.9p . 

Then there is a constant C = C{q) > such that 

(6-4) \Li.A^)\ < ra(T)||fc'+ifl('+i)||^^(jg^^j_,^^, 

where (A:9+iu(9+i))i(t) = kl'/^^Wu^^^'^^^ fort £ 1,^, and where the stability 

factor S{T) is given by S(T) — WJ'^ (tt1^^u,u, dt. Furthermore, if assump- 

tions (Al)-(A5) hold, then there is a constant C — C{q, Ck, a) > such that 

(6.5) |L^,,(e)| < C5(r)||fc2^+i|l(2^+i)||^^(j^^^,^^^^, 

where {k^i+^u'-^i+%(t) = k^-^'^^Cf' \\u'f''~^^''\\L^(u^) for t G hj , and where the 
stability factor S{T) is given by 

S{T)= / C/||$|U^(r,i^)dt = ^if„Q||$|U^(r„,z^). 
Jo . 1 
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Using the stability estimate proved in section 21 we obtain the following bound 
for the stability factor S{T). 

Lemma 6.7. Assume that K^CqCf < 1 for all time slabs Tn, with Cq > the 
constant in Theorem 14.21 and take g = in (|2.9I) . Then 

(6.6) S{T) < ||Vlke^''5,T^ 

where Cf = max[o_T] Cf. 

Proof. By Theorem 14. 2[ we obtain 



M 



\mLUT„,i^) < exp ( J < c< 

and so 

n—1 n—1 

< Mi^ I CqCfe^'-^'fUt < Mi^e^^'^f^. □ 



C,Cf{T-T„) 

loo I 



\m—n+l 



M M 



,C,Cf{T-T„) 



Finally, we rewrite the estimates of Theorems 16.51 and 16.61 for special choices of 
data ■0 and g. We first take V' = 0. With (7„ = for n^i^ gi{t) = for f ^ 7^, and 

9i{t) = sgn{ei{t)) / kij , telij, 

we obtain L^ g{e) — y~ Ii ^^'^ ll^illi=o(-fii) — ^-^ii-gi^) ^"7 an inverse esti- 

mate. By definition, it follows that ||ei||i^(/;^) < CL^^g{e) + Ck1y^^\\uf^~^^\\L^(^j..y 
Note that for this choice of g, we have |15||li([o,t],;2) = \\9\\li{[q.t],i^) = 1- 

We also make the choice 5 = 0. Noting that e{T) = e(T), since nuiT) — u{T), 
we obtain 

L^,,(e) = (e(T),0) = |e,(r)| 

for ipi = sgn(ei(r)) and V'n = for n ^ i, and 

L^,,(e) = (e(r),V) = ||e(r)||,, 

for = s{T) /\\e{T)\\i2. Note that for both choices of if;, we have HV'llioo ^ 1- 

With these choices of data, we obtain the following versions of the a priori error 
estimates. 

Corollary 6.8 (a priori error estimate for mcG{q)). Let U he the mcG(q) 
solution of (jl.ip . Then there is a constant C = C{q) > such that 



(6.7) ||e|U^([o,T].^)<C5(T)||fc^+is('+i)| 



ioo([o,T],;2) 



where the stability factor S{T) = \\J^ (tt'^^qU, u, •)$||;2 dt is taken as the maximum 
over = and \\g\\Li{[o.T],i^) = 1- Furthermore, if assumptions (Al)-(A5) and the 
assumptions of Lemma \6.7\ hold, then there is a constant C — C{q,Ck,a) such that 

i=e([0,T]Jl) 



||e(r)||,^<ra(T)||fc2?s(2?)| 
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for p — 2, oo, where the stability factor S{T) is given by S{T) = e'~^''^f'^ . 

Corollary 6.9 (a priori error estimate for mdG((7)). Let U be the mdG(g) 
solution of Then there is a constant C = C{q) > such that 



(6.9) ||e|U^([o,T].;^) <C5(r)||A:''+iil('+i)| 



ioo([o,T],;2)' 



where the stability factor S{T) ~ || J^(7r|jQU, u, •)<&||;2 dt is taken as the maximum 
over ij: = and \\g\\Li{[Q,T].i^) — 1- Furthermore, if assumptions (Al)-(A5) and the 
assumptions of Lemma \6.7\ hold, then there is a constant C = C{q,Ck,a) such that 



(6.10) l|e(T)||z^ < C5(r)||fc29+is(29+i)| 



L^([O.T],h) 



for p — 2, oo, where the stability factor S{T) is given by S{T) — e*^'*^^"^. 

The stability factor S{T) that appears in the a priori error estimates is obtained 
from the discrete solution $ of the dual problem (|4.ip . and can thus be computed 
by solving the discrete dual problem. Numerical computation of the stability factor 
reveals the exact nature of the problem, in particular, whether or not the problem is 
parabolic; if the stability factor is of unit size and does not grow, then the problem is 
parabohc by definition; see [36] . 

6.3. A note on quadrature errors. The error representations presented in 
section 16.11 are based on the Galerkin orthogonalities of the mcG{q) and mdG(g) 
methods. In particular, for the mcG{q) method, we assume that 



[ {R{U,-),^)dt^O. 
Jo 



In the presence of quadrature errors, this term is nonzero. As a result, we obtain an 
additional term of the form 

^{f{U,-)-f{U,-),^)dt, 

where / is the interpolant of / corresponding the quadrature rule that is used. A 
convenient choice of quadrature for the mcG(g) method is Lobatto quadrature with 
g + 1 nodal points [32], which means that the quadrature error is of order 2{q+l) — 2 = 
2q and so (super)convergence of order 2q is obtained also in the presence of quadrature 
errors. Similarly for the mdG(g) method, we use Radau quadrature with q + I nodal 
points, which means that the quadrature error is of order 2{q + 1) — 1 — 2q + I, and 
so the 2q + I convergence order of mdG((7) is also maintained under quadrature. 

7. A numerical example. We conclude by demonstrating the convergence of 
the multiadaptive methods in the case of a simple test problem. 
Consider the problem 

ill = U2, 

W3 = -1*2 + 2m4, 
Ui — ui — 2i*3, 

1*5 = —1*2 — 2l*4 + 4m6, 
-1*6 = "1 + 21*3 — 4U5 
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Fig. 6. Convergence of the error at final time for the solution of the test problem 117.11 1 with 
mcG(g) and mdG{g), <? < 5. 



Table 1 

Order of convergence p for mcG{q). 



mcG((ji) 


1 


2 


3 


4 


5 


P 


1.99 


3.96 


5.92 


7.82 


9.67 


2q 


2 


4 


6 


8 


10 



Table 2 

Order of convergence p for mdG(g). 



mdG(g) 





1 


2 


3 


4 


5 


P 


0.92 


2.96 


4.94 


6.87 


9.10 




2q + l 


1 


3 


5 


7 


9 


11 



on [0, 1] with initial condition u{0) = (0, 1, 0, 2, 0, 3). The solution is given by u{t) = 
(sint, cost, sint + sm2t, cost + cos 2t, sint + sin2t + sin4t, cost + cos 2t + cos4i). For 
given ko > 0, we take ki{t) = fco for i = 1,2, ki{t) = fco/2 for z = 3, 4, and fci(t) = k^/A 
for i = 5, 6, and study the convergence of the error ||e(T)||/2 with decreasing fco- From 
the results presented in Figure |6] and Tables [1] and [2l it is clear that the predicted 
order of convergence is obtained. 
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